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Abstract 

We consider a model of bosons on a regular lattice with a kinetic energy 
due to hopping among sites and a potential energy due to strong on site in- 
teraction. A superfluid phase is expected when the ground state of the local 
energy is doubly degenerate. We consider a new scheme of simmetry break- 
ing associated to the superfluid phase in which the order parameter is the 
statistical average of the quantum coherence operator associated to the su- 
perposition of the degenerate local ground states. In the strong coupling limit 
a systematic expansion of the free energy can be performed in terms of the 
hopping amplitude at constant order parameter. Within such an expansion 
we obtain a self-consistent equation for the order parameter. The first order 
approximation gives, in the case of degeneracy between single occupied and 
empty state, the same result of the standard mean field approximation for 
the "hard core bosons" . This new approach to the superfluid phase is shown 
to have a natural application to the implementation of quantum computation 
on a superfluid. 
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The achievement of Bose Einstein Condensation experiments and the interest to 

the quantum computing lead naturally to the proposal of using superfluid macroscopic states 
as qubits On the other hand trapping of bosons in optical lattices, P opens the possi- 
bility of new gating procedures. It is therefore of interest to improve the theory of bosonic 
particles with strong repulsive interaction and hopping among sites of a regular lattice and 
to investigate local excitations as possible qubits. 

It is well known that at zero temperature the transition from a superfluid to a Mott 
insulator phase is expected at a critical value of the ratio between the interaction strength 
and hopping amplitude except for integer values of the ratio between the chemical potential 
and the interaction strength = Un* 1^,0]. These values correspond in the atomic limit, i.e. 
for vanishing hopping amplitude, to a degeneracy of the ground state which corresponds to 
the same statistical weight of states involving n* and n* + 1 particles. The associated energy 
can be assumed to vanish for a suitable choice of the zero energy level. A strong coupling 
expansion around the atomic limit has been proposed in ||^. This approach depends however 
on a "factorization" assumption. We show that an alternative approach can be introduced 
which avoids any factorization, and gives, to the first order in the hopping amplitude, the 
classical mean field result. The key point is a new definition of the superfluid phase in which 
the order parameter is associated to the statistical average of suitably defined local quantum 
coherence operator. 

We define the coherence of a system with respect to a particular state as the probability 
of finding the system, or part of it, in that state in the ensemble of all possible states. If 
the system is defined on a lattice, as in our case, and the ensemble is that of the statistical 
equilibrium, it is possible to define a local quantum coherence probability ( LQCP ) as the 
average of the projection operator iVi = \Ai >< Ai\ associated to the local state. 



H is the energy operator in the grand canonical ensemble. It is the sum of a local repulsive 
term Hq = — ^) — fJ'J^i'iT'i { U > ) and a kinetic one Hi = —JJ^ik^ikbfbk 
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(1) 



where J > and eij is equal to if site i is not next neighbour of site j and 1 otherwise. 
We consider here |y4j > as a superposition of the two degenerate groundstates of the local 
energy operator. In the subspace of two degenerate groundstates it is useful to define a Pauli 
vector operator whose components are 

CTji = {\n* + 1 >< n*\ + \n* >< n* + 1|) 
= - (In* + 1 >< n*\ - In* >< n* + 1|) 

= {\n* + I >< n* + l\ - \n* >< n*\) (2) 

The identity restricted to this subspace is defined as 

= (|n* + 1 >< n* + 1| + \n* >< n*\) (3) 

Analogously to the spin case the LQC operator can be expressed in terms of the restricted 
identity and the scalar product of a Pauli operator with the unitary vector { i^i )■ 

TTj = ^CTio + (Tj ■ Z/j (4) 

We show that superfluidity can be associated to a non trivial LQCP ( p{TTi) 7^ | )• Differently 
from the classical approach the gauge symmetry breaking is achieved introducing a field 
defined in the degenerate groundstate subspace in the 1-direction. It is then natural to 
assume this direction as polar axis and thus tTj becomes 
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TTi = - 
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o-jo + cos^jCTji + sin^i (e*'^'(Tj +e '''''af^ (5) 



Raising and lowering operators af = a,;2 izajs appear in (||). Spontaneous gauge symmetry 
breaking is studied in terms of the generating function associated to the 1-direction in the 
degenerate groundstates subspace. 



G = tr 



(6) 



The local order parameter m is defined as mj = ^ In G. Symmetry breakdown occurs if the 
order parameter does not vanish in the limit of Aj — > A — > 0. The next step is to introduce 
the free energy as Legendre transform of In G 



A = \nG-Y.\m, (7) 
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The free energy A depends on rrii instead of the symmetry breaking amphtude Aj, and A 
becomes a function of J determined by the condition 

Note that now the order parameter rrij does not depend on the hopping amphtude J and 
thus can be fixed from the relation with A vahd at J = 

tr {e-f^^exp Ai(0)o-ii] an} 



--<aii> (9) 



tr{e-^^exp[J:^X^iO)au]} 
Following the work of we introduce an expansion in the reduced hopping amplitude (3 J. 
The difficulty is due to the presence, unlike the Ising case, of an atomic Hamiltonian which 
does not commute with the hopping one. Standard formalism of field theory is not suitable 
for the expansion of the generating function of (|^). It is worth to note that if the trace is 
performed on the eigenstates of the atomic Hamiltonian, only the degenerate groundstates 
are affected by the symmetry breaking operator. The expansion in power of (3 J around the 
value (3J = of free energy A can be expressed as follows 

A = Ao + Ai (10) 

Where Aq is defined as 

Ao = In tr |e-^^«exp ^ Ai(0) (a^ - 
It is immediately shown that 

= E In {2 [cosh (A,(0)) - 1^)]] - m,A,(0) (12) 

i 

where '~i{(3U) = 1 — | J2k^^^^''- It is important to note that in the strong coupling limit ( 
f3U ^ 1 ) 'j{l3U) vanishes. In this limit the result valid for the hard core boson model is 
recovered. Ai can be calculated as power expansion in the reduced hopping amphtude /3J 



(11) 
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A^ = A'{PJ) + A"^^ + --- (13) 

where 

A' = J2(^ijmimk{n* + 1) (14) 

ik 

while A" has been calculated up to the first order in 

The local order parameter rrii can be evaluated by the extremum condition for free energy 
= 0. Taking into account the equation (|) for /3J = 

Since rrii is fixed equal to < an > for any reduced hopping amplitude (3 J, it is in particular 
equal to < an > when f3J = 0, which gives us the important relation 

' cosh(Ai(0)) + 7(/3t/) ^ ^ 

From (0) and (|T7D and from the homogeneity hypothesis on the interaction metrix ( e^j = 1 
or ) we obtain a self-consistent equation for the order parameter rrii = ^n. Taking into 
account in Ai only the first order we obtain 

^ ^ sinh(2JD/3m(n* + l)) 

cosh {2JDPm{n* + 1)) + -f{pU) ^ ' 

Where D is the number of next neighbours sites. In the strong coupling limit we can neglect 
7 in (|T8|) and when n* = we obtain the classical result of the mean field approximation 
for the hard core bosons model. As we shall see in the following, the knowledge of the 
expansion up to the second order in the reduced hopping amplitude (3 J becomes important 
to determine the relaxation of local perturbation 0. 

As far as the behaviour of the LQCP is concerned we note that in the non superfiuid 
phase only the statistical average of aio does not vanish and then LQCP becomes trivial 



{pi{ni) = |). On the other hand it is easily seen that in the superfluid phase the sponta- 
neous symmetry breaking imphes a non trivial LQCP pi(7r,) = i (1 + nii cos 9i) . Following 
the work of reference it is convenient, for the implementation of a quantum computer on 
a superfluid phase, to identify the 0-logical state in the i-th site with the dependence of pi 
on cos 9i. To obtain a 1-logical state we need a local perturbation of LQCP. We consider the 
effect of an istantaneous local variation of the chemical potential asssociated to a pertur- 
bation Hamiltonian Hp{t) = 6{t — 0^)Hp with Hp = ^iUiOi. The time evolution of LQCP 
can be obtained in terms of the time evolution of the statistical average of time dependent 
raising and lowering operators. 

~ tr(e-/3^^exp(E.Afeai,)) ^^^^ 

Where cxf (t) are operators in the Heisenberg representation. Simple results are obtained 
assuming that the impulsive perturbation acting on a given site has no influence on next 
neighbours sites, and if the evolution can be approximated to the first order in the reduced 
hopping amphtude. 

lnpf{t) = lnp,^(0+) ± ithi (20) 

where hi = Jy/n* + 1 J2k ''^k^ik has the meaning of the internal static field if m is the static 
order parameter. At this stage of approximation we obtain a free "rotation" around the 
symmetry breaking direction completely analogous to that found in the magnetic case. We 
expect that decoherence effects will appear when the correlation of different sites degrees of 
freadom are taken into account. An evaluation of decoherence can be obtained introducing, 
in analogy with the reference 0, a dynamical function 

Gf{X,t)=\ni^Tr e-^"" exp (^{±tJt^/^^^F+le,k + Xk)(rk?j | (21) 

We assume that pf{t) can be written as 

pf{t)= sin 9iAf{t) (22) 
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The dependence on sin 61 can be considered as the manifestation of the presence of a 1-logical 
states in site i 0. The time behaviour of 1-logical states Af{t) are given in terms of the 
generating dynamical function. 

Af{t) = exp {Gf{t) - Gf{0)) ^ (23) 

The equation ( p3D is actually an extrapolation of the result obtained in the first order expan- 
sion in the reduced hopping amplitude {(3 J). The calculation of the dynamical generating 
function can be done introducing a Legendre transform at fixed time and reduced interaction 

Ftim, t) = Gf{X, P,t)-Y^ Afc(/9, t)mk (24) 

k 

Following an analysis similar to that of reference we obtain 

Ft{m,(3,t) ^ Fi{m,(3)±tJtY, ei^nik - Tf (25) 

k 

where T is equal to J^(n* + iyj2k^ik{^ — '^^) + 0{f3U). We see that within the present 
approximation decoherence decreases as the superfiuid parameter saturates (m^ 1). 

In conclusion we have studied the superfiuid phase of a boson system on a regular lat- 
tice with a strong on site interaction. The degenerate groundstates of the energy operator 
associated to each site are used to define a qubit. Hopping of particles among sites defines 
the interaction of the qubits with the physical environment. The superfiuid phase is defined 
in terms of a generating function associated to the local coherence operator. This operator 
defines the direction in the subspace of the degenerate groundstates along which symmetry 
breaking occurs. Non trivial properties of LQCP are found in the superfiuid phase. Quan- 
tum gating on a single qubit is associated to the modification of LQCP due to a suitable 
local impulsive gauge perturbation. Decoherence appears only when correction to the mean 
field approximation are taken into account. Slowing down of decoherence is obtained as the 
superfiuid order parameter approaches saturation. In other words this particular environ- 
ment gives rise, in the superfiuid phase, to an effective field along which qubits tend to align 
( 0- logical states ). Excitations of qubits in the orthogonal direction ( 1-logical states ) will 
relax slower and slower as superfiuid order parameter saturates. 
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